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I The multiple disorder problem seeks to determine a sequence of stopping times which 

, are as close as possible to the unknown times of disorders at which the observation process 

I changes its probability characteristics. We derive closed form solutions in two formula- 

tions of the multiple disorder problem for an observable Brownian motion with switching 
constant drift rates. The method of proof is based on the reduction of the initial problems 
to appropriate optimal switching problems and the analysis of the associated coupled free- 
boundary problems. We also describe the sequential switching multiple disorder detection 
procedures resulting from these formulations. 
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H ! 1. Introduction 

Suppose that at time t = we begin to observe a sample path of some continuous process 
X = {Xt)t>Q with probability characteristics changing at some unknown disorder times (?7„)neN 
at which an unobservable two-state process O = {Qt)t>o switches between one state and the 
other. The switching multiple disorder problem is to decide at which time instants ('r„)„gN 
one should give alarm signals to indicate the occurrence of changes in the current state of the 
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process 9 as close as possible to the initial disorder times (?7n)neN- Such disorder (or change- 
point) detection problems have originally arisen and still play a prominent role in quality 
control, where one observes the output of a production line and wishes to detect deviations 
from acceptable levels. After the introduction of the original control charts by Shewhart [29] , 
various modifications of the disorder problem have been recognized (see, e.g. Pages [23]) and 
implemented in a number of applied sciences (see, e.g. Carlstein, Miiller and Siegmund [T2]). 

The problem of detecting a single change in the constant drift rate of a Brownian motion 
(Wiener process) was formulated and explicitly solved by Shiryaev [30]- [31] and [31]- [35] (see 
also Shiryaev [361 Chapter IV] and Peskir and Shiryaev [2H Chapter VI, Section 22] for further 
references). The optimal time of alarm was sought as a stopping time minimising a linear 
combination of the false alarm probability and the average time delay in detecting of the disorder 
correctly. Shiryaev pO] and [32] also proposed another formulation of the problem in which the 
occurrence of a single change should be preceded by a long period of observations under which 
a stationary regime has been established. The resulting optimal multistage detection procedure 
consisted in searching for a sequence of stopping times minimising the average time delay given 
that the mean time between two false alarms is fixed. More recently, Feinberg and Shiryaev 
[TB] derived an explicit solution of the quickest detection problem in the generalized Bayesian 
formulation and proved the asymptotic optimality of the associated detection procedure for the 
related minimax formulation. Extensive overviews of these and other related sequential quickest 
change-point detection methods were provided in Shiryaev [37] and Poor and Hadjiliadis [26] . 

In the present paper, we formulate and solve the switching multiple disorder problem for 
an observed Wiener process X changing its drift rate from fij to /ii-j when O changes its 
state from j to 1 — j , for every j = 0, 1 . In contrast to the problem of detecting a single 
change, in the switching multiple disorder problem, one looks for an infinite sequence of the 
alarm times (T„)„gN minimising a series of linear combinations of discounted average losses due 
to false alarms and delay penalties in detecting of the disorder times (77„)„eN correctly. We 
propose two different formulations of the problem based on a specification of dynamics of the 
process G. In the first formulation, B is assumed to be a continuous time Markov chain of 
intensity A, the dynamics of which are not influenced by the alarm times (r„)„gN- In the second 
formulation, it is assumed that the subsequent time rjn, at which B changes its state, can only 
occur after the previous alarm is sounded at r„„i. Moreover, it is assumed that the differences 
iVn — Tn-i)neN form a sequence of (conditionally) independent exponential random variables. 
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Apart from other possible areas of application, such a situation usually happens in models 
of liquid financial markets having trading investors of different kinds. It is natural to assume 
that the small investors can only influence little fluctuations of the market prices of risky assets, 
while the large investors can affect the pricing trends as well, by means of either buying or selling 
substantial amounts of assets. More precisely, the pricing trends should either rise up or fall 
down at some random times, after essential amounts of assets are bought or sold, respectively. 
We can thus consider a model of a financial market of such kind in which the dynamics (of the 
logarithms) of the asset prices are described by a Brownian motion with switching drift rates. 
We may further assume that our model allows for an infinite number of transactions (free on 
charge) on the infinite time interval and use an exponential constant discounting rate r , which 
can be chosen equal to the riskless short rate of a bank account. The problem of detecting of a 
single change in the probability characteristics of accessible financial data, which is associated 
with the appearance of arbitrage opportunities in the market, was considered by Shiryaev P?] . 

In the present paper, we reduce the initial multiple disorder problems to appropriate optimal 
switching problems for filtering estimates of the current state of the unobservable drift rate of 
a Brownian motion. The use of exponential discounting makes our problem well connected to 
the problem of single disorder detection with exponential delay penalty costs studied by Poor 
|25j . Beibel [S], and Bayraktar and Dayanik [3J. We show that the optimal switching times 
can be expressed as the first times at which the appropriate posterior probability processes hit 
certain constant boundaries. We derive closed form expressions for the resulting Bayesian risk 
functions and the optimal switching boundaries by means of solving the associated coupled free- 
boundary problems for ordinary differential operators. We also construct sequential switching 
multiple disorder detection procedures resulting from the two formulations. 

Optimal switching problems represent extensions of stopping problems and games in which 
one looks for an infinite sequence of optimal stopping times. A general approach for studying 
such problems was developed in Bensoussan and Friedman [S]-[in] and Friedman [T7] (see also 
Friedman [TSl Chapter XVI]). This investigation was continued by Brekke and 0ksendal [TT] , 
Duckworth and Zervos ^4J , Yushkevich and Gordienko [IQ] , and Hamadene and Jeanblanc [19] 
among others for the continuous time case, and by Yushkevich [38]-[39| for the discrete time 
case. Other optimal switching and impulse control problems involving hidden Markov chains 
in the observable jump processes were recently studied by Bayraktar and Ludkovski [S]-IZ]- 

The paper is organized as follows. In Section 2, for the initial multiple disorder problems, we 
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construct the appropriate optimal switching problems and reduce the latter to their equivalent 
coupled optimal stopping problems. In Section 3, we derive closed form solutions of the associ- 
ated coupled free-boundary problems, which are expressed in terms of Heun's double confluent 
functions and Kummer's confluent hypergeometric functions. In Section 4, we verify that the 
solutions of the coupled free-boundary problems provide the solutions of the initial optimal 
switching problems, and describe the resulting sequential switching multiple disorder detection 
procedures. The main results of the paper are stated in Theorems 4.1 and 4.2. The optimal 
sequential detecting schemes are displayed more explicitly in Remark 4.3. 

2. Formulation of the problems 

In this section, we give two formulations of the switching multiple disorder problem for an 
observed Brownian motion (see, e.g. [361 Chapter IV, Section 4] or [211 Chapter VI, Section 22] 
for the single disorder case). In these formulations, it is assumed that one observes a sample 
path of the Brownian motion X with the drift rate switching between /iq and fii at some 
random times {rin)neN- 

2.1. (The setting.) Let us assume that all the considerations take place on a probability 
space (f2, Q, Pt^) with a standard Brownian motion (Wiener process) B = {Bt)t>o started at 
zero under P^^ • Suppose that there exists a right-continuous process G with two states and 
1, having the initial distribution {1 — 7r,7r} under P^^, for tt G [0,1]. It is assumed that the 
process G is unobservable, so that the switching times = inf{t > rj^-i \ Qt 7^ 0t?„-i}) fo^' 
n G N, with r^o = 0, at which G changes its state from j to 1 — j, for every j = 0, 1, are 
unknown, that is, they cannot be observed directly. 

Suppose that we observe a continuous process X = {Xt)t>o solving the stochastic differential 
equation: 

dXt= {i2o + {121-120) et)dt + a dBt (Xo = 0) (2.1) 

where /xq 7^ /^i and a > are some given constants. Being based upon the continuous 
observation of X , our task is to find among non-decreasing sequences of stopping times (r„)„gN 
of X (i.e., stopping times with respect to the natural filtration J^t = cr{Xs | < s < t) of the 
process X, for t > 0) at which the alarms should be sounded as close as possible to the 
unknown switching times of the process G. More precisely, the switching multiple disorder 
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problem consists of computing the Bayesian risk functions: 



oo 

inf 



a e 



-'"ro,2fc- 



1 /(e 



''0,2*; -1 



be 



0) 



(2.2) 



F;(7r)= inf 



/-fO,2fc-l+J 

+ E/ e-''l{et=j)dt 



5e-rn,2fc_i /(e,^ ^^_^ = 0) + ae-"^i-2'= 7(9 



Tl,2fc 



(2.3) 



rTl,2k-l+j 

+ E / e-^*J(e, = l-j)c^t 

j=0 "'n,2fc-2+i 

and finding the non-decreasing sequences of optimal stopping times (T*„)„gN such that r*g = 0, 
i = 0, 1, at which the infima in fl2.2p and (12.31) are attained, respectively, where /(■) denotes 
the indicator function. Note that the function V*{tt) expresses the Bayesian risk of the whole 
sequence (rj„)„gN in the case in which the process 6 starts at Oq = 1 — "i, for every i = 0,1 
fixed, and all tt G [0, 1] . We therefore see that i?^ [e"*"^'-" /(O^-^ ^ = j)] is the average discounted 
loss due to a false alarm, and E^, [ JJ^"'" ^ e"*"* /(©* = 1 — j) dtj expresses the average discounted 
loss due to a delay in detecting of the time at which G changes its state from j to 1 — j 
correctly, for every i,j = 0, 1 and any n G N. In this case, a,b > are costs of false alarms 
and r > is a discounting rate. 

Using the fact that (rj_„)„gN is a non-decreasing sequence of stopping times with respect to 
the filtration {J^t)t>o, by means of standard arguments, which are similar to those presented in 
[361 pages 195-197], we get that: 

[e— - J(e.,,„ = j)] = E^ [E^ [e— /(e.,,„ = j) \ T^J] = E^ [e— - P.(e.,,„ = j \ J".,, 
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holds for every z, j = 0, 1 and any n eN. 

We further consider two different formulations of the problem, depending on the specified 
dynamics of the process G. The first formulation does not involve any influence of the alarm 
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times Ti^n on the times of changes //„ . In the second formulation, it is assumed that the change 
at rin can occur only after the previous alarm is sounded at Tj „_i, for every i = 0,1 and any 
n e N. 

2.2. (The first formulation.) Suppose that G is a continuous time Markov chain which 
is independent of the Brownian motion B and has the initial distribution {1 — 7r,7r} under 



P^. Assume that G has the transition-probability matrix {e 



-At 



} and 



the intensity-matrix {—A, A; A, —A}, for all t > and some A > fixed. In other words, the 
Markov chain G changes its state at exponentially distributed times of intensity A , which are 
independent of the dynamics of the Brownian motion B . Such a process G is called telegraphic 
signal of intensity A in the literature (see, e.g. [211 Chapter IX, Section 4] or (TS] Chapter VIII]). 

It thus follows from [21] Chapter IX, Theorem 9.1] (see also [21; Chapter IX, Example 3]) 
that the posterior probability process 11 = (n()(>o defined by Ut = P-K^Qt = 1 1 J^t) solves the 
stochastic differential equation: 



a 



dUt = A(l - 2Ut) dt 
where the innovation process B = (Bt)t>o defined by: 



Util-Ut)dBt (no = 7r) 



(2.6) 



Bt = ^ (^Xt - J (yUo + (/Ui - /io) n^) ds^ 



(2.7) 



is a standard Brownian motion according to P. Levy's characterization theorem (see, e.g. 
Chapter IV, Theorem 4.1]). It is also seen from (12. 6p that 11 is a (time- homogeneous strong) 
Markov process with respect to its natural filtration, which obviously coincides with {J^t)t>o- 

Taking into account the expressions in (12.41) and (12. Sp . we therefore conclude that the 
Bayesian risk functions from (12. 2 p and (12. 3p admit the representations: 
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where the infima are taken over all sequences of stopping times (Tj^„)„gN, ^ = 0, 1 , of the process 
n. By virtue of the strong Markov property of the process 11, we can reduce the system of 
(12. 8p and (12. 9 p to the following coupled optimal stopping problem: 



ae-^^° + / e-'^* (1 - Ut) dt + V*{UJ 
Jo 



'-T1 I 



(2.10) 
(2.11) 
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where the infima are taken over all stopping times r,, i = 0, 1, of the process 11 with P7r(no = 
tt) = 1 . We further search for optimal stopping times in (I2.10p and (12. lip of the form: 

T* = M{t > I < and t* = M{t >0\Ilt>K} (2.12) 

for some < < /i* < 1, where g^, is the largest and /i* is the smallest number vr from [0, 1] 
such that V'o*(7r) = an + V\*(7r) and V{{tt) = 6(1 - tt) + V'o*(7r) holds, respectively. This fact 
implies that the sequences of stopping times (T*„)„gN given by: 

T*2k-i+i = inf{t > r*2fc_2+i I < g^} and r*2;._. = inf{t > rl^^-i-i I > K} (2.13) 

for every i = 0,1 and any A; G N, are optimal in the problems of (12. 8 p and (12. 9p . 

2.3. (The second formulation.) As that is the case in the previous formulation, for every 
i = 0, 1, let us denote by (Cj,2fc-i)fceN and {Ci,2k-i+i) ken the sequences of alarm times sounded 
to indicate that the state of G has been changed from to 1 or from 1 to , respectively. Let 
us now assume that the switching time rjn of the process G can only occur after the previous 
alarm is sounded at O.n-i, for any n G N. Suppose that (^j,n)nGN defined by = Vn — Ti,n-i 
forms a sequence of (conditionally) independent non-negative random variables such that ^i^n is 
independent of the Brownian motion B on the time interval [Cj,n~i, Ci,n] ■ Moreover, we assume 
that the properties Pj^iVn = Ci,n-i | ^o,n-i) = ^Q,n-i and Pn{r]n > t\r]n > Ci,n-i, ^o,n-i) = 
g-A{t-Ci,n-i) hold for all t > Ci,n-i and some A > fixed. In other words, the process G 
changes its state in the exponential time C,i,n = Vn ~ Ci,n-^i of intensity A after the time of the 
previous alarm O.n-i, where ^i^n does not depend on the subsequent fluctuations of the process 
B. 

It thus follows from pT; Chapter IX, Theorem 9.1] (see also [2T1 Chapter IX, Example 2] or 
[151 Chapter VIII]) that the posterior probabihty process 11 solves the stochastic differential 
equation: 

duf'^ = -xui'^dt + ^^^^uf\i-uf^)dB, (n(°) =n?) ) (2.14) 



for Ci,2k-2+i <t < Ci.,2k-i+i and 
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(2.15) 



foi' Ci,2fc~i-i l£ t < Ci,2fc-j5 where the process _B is defined in (12.71) and turns out to be a standard 
Brownian motion on the time intervals [0,n-i; Ci,n] , for every z = 0, 1 and any fc, n G N. 

Taking into account the expressions in (12 ■4p and (12. 5p . we may conclude that the Bayesian 
risk functions in this formulation are given by: 
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where the infima are taken over all sequences of stopping times (Cj,n)ngN of the processes 
n(*) = (nf^)t>o, i = 0,1, solving the stochastic differential equations in (12.141) and (12.151) . 
respectively. By virtue of the strong Markov property of the processes n», i = 0,1, we can 
reduce the system of (12.161) and (12.171) to the following coupled optimal stopping problem: 

r-Co 
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where the infima are taken over all stopping times Q of the processes H^*) , i = 0,1, respectively. 
We further search for the optimal stopping times in (12.181) and (I2.19P of the form: 



Co = inf{t > I < and Ci* = mf{t > | U\'> > 



(1) 
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for some < < < 1, where is the largest and q^: is the smallest number vr from [0, 1] 
such that Uq{7t) = an + U*{tt) and f/i (vr) = 6(1 - tt) + Uq{h) holds, respectively. This fact 
implies that the sequences of stopping times (C*n)neN given by: 

C- = inf{t > C- 2fc-2+. I < P*} and C* 2.-. = inf{t > C* 2^-1-^ I ^i'^ > Q*} (2-21) 
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for every i = 0,1 and any k eN, are optimal in the problems of f l2.16p and fl2.17p . 

Remark 2.1. Recall that (Cj,n)ngN is a non- decreasing sequence of stopping times with 
respect to the filtration {J^t)t>o- Then, by virtue of the assumption that rjn > Ci,n-i holds, 
standard arguments show that the equalities: 



e-'' I{Qt = j) dt 



Ci , n — 1 



= Et, 
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e /(Ci,n-1 <t,r]n<t,t< 0,n) dt 
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dt 



are satisfied for every i,j = 0,1 and any n G N. This fact builds a connection between 
the introduction of exponential discounting into the switching multiple disorder problem of 
(I2.16p -( l2.17p and the consideration of single disorder detection problems with exponential delay 
penalty costs studied in [25], [8] and [3]. 

2.4. (Coupled free-boundary problems.) Standard arguments based on an application of 
Ito's formula (see, e.g. [2Q- Chapter V, Section 5.1] or [22; Chapter VII, Section 7.3]) imply that 
the infinitesimal operator L of the process 11 from (12.60 acts on an arbitrary twice continuously 
differentiable (locally) bounded function E{tt) according to the rule: 

(LF)(7r) = A (1 - 2n) F'{n) + 1 - F"{n) (2.23) 

for all vr G (0, 1). In order to find the unknown value functions V"Q*(7r) and V*{tt) from fl2.18p 
and (I2.19P as well as the unknown boundaries g^, and /i* from (12.120 . we may use the results 
of the general theory of optimal stopping problems for continuous time Markov processes (see, 
e.g. [361 Chapter III, Section 8] and [211 Chapter IV, Section 8]). More precisely, we formulate 
the associated coupled free-boundary problem: 
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Vo(vr) = an + Vi{n] 


for n < g, 


V,in) 


= b{l- 
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for n > h 


(2.27) 




Vo{7r) < an + Vi{n] 


for n > g, 


V^in) 


< 6(1 - 


vr) + K)(vr) 


for n < h 


(2.28) 


(LFo 


-rVo)in)>-{l-n] 


for n < g, 


(L^i 


-rVi){n 


) > —n for 


n > h 


(2.29) 



with < g < h < 1, where the instantaneous-stopping and smooth-fit conditions of fl2.25p and 
fl2.26p are satisfied at g^ and h^. Note that the superharmonic characterisation of the value 
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function (see, e.g. p6j Chapter III, Section 8] and [Ml Chapter IV, Section 9]) implies that 
Vq{tt) from (12.101) and V{{tt) from (12. lip are the largest functions satisfying the expressions 
in fl2:2l|) - (12:251) and (12:271) - (12:281) with the boundaries and K. 

Furthermore, standard arguments show that the infinitesimal operator Lj of the process 
n(^) from ( I27[4l) -( 12351) acts on an arbitrary twice continuously differentiable (locally) bounded 
function -F(7r) according to the rule: 

(LoF)(7r) = -XnF'in) + \ - nfF"{7r) (2.30) 

iUF){n) = A (1 - vr) F'in) + ^ - -nf F"in) (2.31) 

for all TT G (0, 1) and every i = 0,1. In order to find the unknown value functions Uq{tt) and 
UKt:) from (I2.18P and (12.190 as well as the unknown boundaries and g* from ( 12.120 . we 
formulate the associated coupled free-boundary problem: 



rf/o)(7r) = -(l- 


vr) 


for TT > p, 


(Lit/i 


— rUi){TT) = —TT for TT < q 


(2.32) 




= a 


p + Ui{p+), 


Uiiq- 


) = h{l-q) + U,{q-) 


(2.33) 




) = 


a + U[{p+), 


U[{q- 


) = -h + U',{q-) 


(2.34) 


Uo{n) = a7T + Ui{ 




for n < p, 


Ui{tt) 


= h{l - tt) + Uq{tt) for TT > q 


(2.35) 


Uo{n) < avr + f/i( 


» 


for TT > p, 


Ui{tt) 


< h {I - tt) + Uo{tt) for tt < q 


(2.36) 


rf/o)(7r)>-(l- 


vr) 


for TT < p, 




— rUi){TT) > —TT for TT > q 


(2.37) 



with < p < q < 1 , where the instantaneous- stopping and smooth-fit conditions of (12.330 and 
(I2.34P are satisfied at and g*. The superharmonic characterisation of the value function 
implies that Uq{tt) from (12.180 and U*{tt) from (I2.19P are the largest functions satisfying the 
expressions in (I2.32p - (l2.33p and (12. 350 - 02. 360 with the boundaries p^ and g*. 



3. Solutions of the coupled free-boundary problems 

In this section we solve the systems of (I2.24p - (l2.29p and (I2.32p - (l2.37p and prove the exis- 
tence and uniqueness of solutions of those coupled free-boundary problems associated to the 
corresponding formulations of the switching multiple disorder problem. 

3.1. (Existence in the first formulation.) The general solutions of the second order ordinary 
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differential equations in (12 .24^ are given by: 

where Cjj, j = 0, 1, are some arbitrary constants, and the functions (5i(vr), i = 0,1, are given 
by: 

Q,(7r) = v/^KT^expf-^ + ii^^) i7Y(-l)'+V,V^,0,e;-^') (3.2) 

Vp(l-7r) pvr y V 1-271 J 

for all vr € (0, 1) with 

P=(— — , ^ = —, V' = ^ + </'- — -1 and ^ = 4:^p-ip. (3.3) 

V a y p 4 p 

Here, the functions Hi{a, f3,'j,6;x) , i = 0,1, are two positive fundamental solutions (i.e. non- 
trivial linearly independent particular solutions) of Heun's double confluent ordinary differential 
equation: 

^ I ~ -Ax^ + 2x + a ^ , + {2a + i)x + S 

Hix)^ -. — r- —: -H (X) H 7 — r- H{x)=0 (3.4) 

^ ' (X-I)3(x+1)3 ^ ^ {X-1Y{X+If ^ ' ^ ' 

with the boundary conditions H{0) = 1 and H'{0) = 0. Note that the series expansion of the 
solution of the equation in (13. 4p converges under all — 1 < x < 1, and the appropriate analytic 
continuation can be obtained through the identity H{a, (3, 7, 6; x) = H{—a, —6, —7, —f3; 1/x) . 
The (irregular) singularities at —1 and 1 of the equation in ( 13. 4p are of unit rank and can 
be transformed into that of a confluent hypergeometric equation (see, e.g. [13] and [28] for 
an extensive overview and further details). According to the results from [271 Chapter V, 
Section 50], we can specify the positive (strictly) convex functions (5i(vr), i = 0, 1, as (strictly) 
decreasing and increasing on the interval (0, 1) and having singularities at and 1 , respectively. 

Taking into account the fact that the function Vo(vr) should be bounded as vr f 1 while 
the function Vi^n) should be bounded at vr J, 0, we must put Cqi = Cio = in (13. ip . Then, 
applying the instantaneous-stopping and smooth-fit conditions from (12.251) and (I2.26P to the 
function in (13.11) . we get that the equalities: 

CuQi{9)-CooQoi9) = M9) and CnQi{h) - CooQo{h) = Ri{h) (3.5) 
CuQ[ig)-CooQ'oi9) = Roi9) and CnQ'M) - CooQW = R[{h) (3.6) 

hold for some < g < h < 1 , where we set: 

1 — 2vr 1 — 2vr 

Ro{n) = -a7r + —— and i?i(vr) = 6 (1 - vr) + (3.7) 
2A + r 2a + r 
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for all IT G [0, 1]. Solving the left-hand part of the system in f l3.5p -( |3^ . we obtain: 

r ( Roig)Q[{9) - K{g)Qiig) . r ( \ - M9)Qoi9) - Ki9)Qo{9) 

Qi{9)Q'oi9) - Q[{9)Qo{9) ''^^^ Qi{9)Q'oi9) - Q[{9)Qo{9) ^ ^ ^ 

and the solution of the right-hand part there gives: 

P. , RmQ\{h) - R[{h)Q,{h) ~ R,{h)Q',{h) - R[{h)Q,{h) 

^ QmQ',{h) - Q',{h)Q,{h) ''^ ^ Qi{h)Q',{h) - Q[{h)Qo{h) 

so that the system in (I3.5l) - (l3.6p is equivalent to: 

Cooi9) = Cooih) and d^iig) = Ci^ih) (3.10) 

for < g < h < 1. It thus follows that the functions: 

Vo{ir;9) = Coo{9)Qo{7^) + ^^7^^^-^ and V,{7r; h) = Cu{h) Q,{n) + (3.11) 

r(/A + r) r{2X + r) 

provide a solution of the system in fl2.24l) - fl2.26l) for any < g < h < 1 fixed. 

3.2. (Uniqueness in the first formulation.) Let us now show that the system in fl3.10p 
with ( I3.8p - (]3.9p admits a unique solution g^ and /;.*. For this, using the standard comparison 
arguments for solutions of the second order ordinary differential equations in (12.240 . we conclude 
that the resulting curves vr i— )• Vo{n; g) and vr i-> Vi(7r; /i) from (I3.1ip do not intersect each other 
on the intervals [g, 1) and (0, h] , respectively, for different < g < h < 1 fixed. We also observe 
by virtue of the properties of the functions Qiiir), z = 0, 1, in (13.20 that Vo(7r; (7) and Vi{7C] h) 
are bounded and concave on [g, 1) and (0, h] , respectively, and such that Vg'(7r; (7) —t- 00 as 
TT I and V{{iT] h) — )■ —00 as vr f 1. On the other hand, using the conditions in (I2.27p . we 
obtain by means of straightforward computations that the inequalities in (12.290 are satisfied 
whenever < g <g and h < h < 1, where we set: 

, and S=l±#4±4 (3.12) 



" 2 + a(2A + r) 2 + 6(2A + r) 

and note that <g < 1/2 < h < 1 holds. Hence, we may conclude that if the conditions: 

V^{g+;g)<a + V;Cg+;h) and V;(h~;h)>-b + V^(h-,g) (3.13) 

are satisfied, then the boundaries g^ and belong to the intervals (0, g) and {h, 1) , respec- 
tively. In other words, the assumptions in (I3.13P describe the set of all admissible parameters 
a,b > for which the free-boundary problem of ( ]2.24p -f l2l29l) admits a unique solution, so that 
the optimal stopping and switching times are given by (I2.12p and (I2.13p . respectively. 
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3.3. (Existence in the second formulation.) The general solutions of the second order 
ordinary differential equations in (12.321) have the form: 

U,in) = Ao + Ai ^^(Tr) + — ^ + + (3.14) 

r(A + r) A + r A + r 

where Dij are some arbitrary constants and the functions Gij{TT), i,j = 0,l, are given by: 

Goo(7r) = (1 - tt) {j^y ^ (7+ - 1, 7+ - 7- + 1; ■^^^) (3-15) 

Goi(vr) = (1 - vr) (y^)'' $ (7+ - 1,7+ - 7- + 1; ;^(^) (3-16) 



and 



GM = vr ( — $ (7+ - 1,7+ - 7- + 1; ^^^^ ) (3-17) 
^11(71) = TT ( ^ (7+ - 1,7+ - 7- + 1; ^^^^^ ) (3-18) 



with 




P = ^^^^ and 7+ = ^ + - ± W U + - + - (3-19) 



cr 



for all 71 G (0, 1). Here, we denote by: 

$(a,/3;x) = 1 + 5:^1^ (3.20) 



IT f $(a,/3;x) $(1 + a - /3, 2 - /3; x) 



^(a,/5;x) (^r(l + a'-'/3)r(/3) r(a)r(2-/3)' J ^^'^^^ 

Kummer's confluent hypergeometric functions of the first and second kind, respectively, for 
/3 ^ 0, -1, -2, . . . and (/3)fc = /3(/3 + l) ■ ■ ■ {(3 + k-l), k eN, where the series in IK20\i converges 
under all x > (see, e.g. [H Chapter XIII] and [H Chapter VI]), and F denotes Euler's Gamma 
function. According to the results from |271 Chapter V, Section 50], we can specify the positive 
(strictly) convex functions Gioi^n), i = 0,1, and Gh^tc), i = 0,1, as (strictly) decreasing and 
increasing on the interval (0, 1) with singularities at and 1, respectively. 

Taking into account the fact that the function Uq{'k) should be bounded as tt f 1 while the 
function f/i(vr) should be bounded at tt J, 0, we must put Ai = -Dio = in fl3.14p . Then, 
applying the instantaneous-stopping and smooth-fit conditions from f l2.33p and fl2.34p to the 
function in fl3.14p . we get that the equalities: 

DnGn{p)-DooGoo{p) = So{p) and Ai G'n(g) - Ao G'oo(g) = 5i(g) (3.22) 
DuG[,ip)-DooGM = S'oip) and DuG[,iq) - DooG'M = S[{q) (3.23) 
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hold for some < p < q < 1, where we set: 

1 — 27r 1 — 27r 
5o(7r) = -a7r + and 5i(7r) = 6 (1 - tt) + (3.24) 

for all TT G [0, 1]. Solving the left-hand part of the system in fl3.22p - fl3.23p . we obtain: 

f)^^^^-^ = 'S'o(p)<^ii(p) - So{p)Guip) Dii{p) - '^o*^^)'^oo(p) - S'q{p)Goq{p) 



G,,{p)G'M - G',,{p)GM ^^^"^ G,,{p)G'M - G[,ip)Gooip) 

(3.25) 

and the solution of the right-hand part there gives: 

. , S,iq)G[,iq) - S[iq)G^r{q) ~ S^iq)G'M ' S[{q)Gu{q) , . 

Gn(g)G"oo(g)-G"n(g)Goo(g) "^^^ G.MG'M ~ G'.MGM ^ ' ^ 

so that the system in fl3.22p - fl3.23p is equivalent to: 

5oo(p) = 5oo(g) and 5n(p) = 5n(g) (3.27) 

for 0<j9<g<l. It thus follows that the functions: 

f/o(7r;p) =Doo(p)G'oo(vr) + ,.\ . ' and [/^(Tr; g) = Z}n(g) G'n(7r) + -— — (3.28) 

r(A -|- rj r(A -|- rj 

provide a solution of the system in f ]2.32p -f !2.34p for any Q < p < q < 1 fixed. 

3.4. (Uniqueness in the second formulation.) Let us finally follow the schema of arguments 
above, to prove that the system of equations in fl3.27p with fl3.25p - fl3.26p admits a unique 
solution and g*. For this, we use the standard comparison arguments for solutions of the 
second order ordinary differential equations in f l2.32p to conclude that the curves vr i— t- Uq{71]p) 
and TT I— 7- f/i(7r;g) from fl3.28p do not intersect each other on the intervals [p, 1) and (0,g], 
respectively, for different < p < g < 1 fixed. We also observe by virtue of the properties 
of the functions Ga^n), i = 0, 1, in f l3.15p and f l3.18p that Uo{n;p) and Ui{n;q) are bounded 
and concave on [p,l) and (0,g], respectively, and such that Uq{tt;p) — )■ oo as vr J, and 
U[{TT;q) —7- — oo as TT 'I" 1. Moreover, using the conditions in f l2.35p . we obtain by means of 
straightforward computations that the inequalities in fl2.37p are equivalent to: 

(2 + a(A + r))p - < -XD^,{q)G[,{p) (3.29) 

X + r 

(2 + b{X + r))(l - g) - < XDoo{p)GUq) (3-30) 

for < p < g < 1. Note that since the derivative G'uln) is positive and increasing from zero to 
infinity, while the derivative G"QQ(7r) is negative and increasing from minus infinity to zero, it is 
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shown by means of standard arguments that the inequahties in fl3.29p and fl3.30p hold whenever 
< p <p and g < g < 1 , where the numbers p and q are set by: 

r 1 , ^ \ + (\ + r)(l + b(\ + r)) 1 

p = pA- TT < - and g = g V — > -• 3.31 

^ ^ (A + r)(2 + a(A + r)) 2 ^ ^ (A + r)(2 + 6(A + r)) 2 ^ ^ 

Here, the couple p and g is determined as a unique solution of the corresponding equations 
instead of the inequalities in fl3.29p and fl3.30p whenever it exists, and p = g = 1/2 otherwise. 
Hence, we may conclude that if the conditions: 

U;,{p+;p) <a + U[ip+;q) and f/((g-; g) > -& + f/o(g-; (3.32) 

hold, then the system in fl3.27p admits a unique solution and g* such that < < p and 
g < g* < 1 . Therefore, the assumptions in f l3.32p describe the set of all admissible parameters 
a,b > for which the free-boundary problem of fl2.32p - fl2.37p admits a unique solution, so that 
the optimal stopping and switching times are given by fl2.2Up and fl2.2ip . respectively. 



4. Main results and proofs 

Taking into account the facts proved above, we are now ready to formulate and prove the 
main assertions of the paper. 

Theorem 4.1. Assume that the conditions in li^.l^) are satisfied with g and h defined in 
lis. 12^) . Then, in the switching multiple disorder problem of Ii2.8\) - I[27^) and Ii2.10\) - f2.11\) for 
the process X from h2. the Bayesian risk functions V*{tt) , i = 0,1, take the form: 

, Vo(n; g^), if g^ < n < 1 

V*in) =^ ' - (4.1) 

air + Vi{n;h^), if < n < g^ 

V*{^)={ (4.2) 
U(l -7r) + 1/o(7r;^,), ifK<TT<l 

and the optimal switching times (r*„)„gN, i = 0,1, have the structure of i\2.13\) . Here, the 
functions Vo{tt; g) and Vi{tt; h) are given by i\3.11\) . and the optimal stopping boundaries g.^ and 
h^: , such that 0<g^:<^<l/2<h<h^:<l, are uniquely determined by the coupled system of 
the equations in 1^3. 10\} with Cii{g) and Cii{h) given by ^3.3i) - ^379i) . where the functions Qiiir) 
and Rii^n) , i = 0,1 , are defined in h3.2\) and \3. 7| ), respectively. 
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Theorem 4.2. Assume that the conditions in ^3. 32\) are satisfied with p and q defined by 
lis. 31\) . where p and q is a unique solution of the system of equations replacing the inequalities 
in Ii3. 29\) - f3. 3(^) whenever it exists, and p = ^= 1/2 otherwise. Then, in the switching multiple 
disorder problem of ^2.16\) - (2.n\ ) and Ii2.18\) - (2.19\) for the process X from ^2. the Bayesian 



risk functions U*{tt), i = 0, 1, take the form: 

jf/o(7r;p*), z/p. <7r<l 

t/o(vr)=< (4.3) 
I avr + [/i(7r; g*), if < tv < p^, 

Ui{'K;q,), < TT < 

U*in) = I (4.4) 

U(l - vr) + f/o(vr;p*), if q* < n < 1 

and the optimal switching times (C*n)neN; i = 0,1, have the structure of Ii2.21\) . Here, the 
functions f/o(7r;p) and f/i(7r;g) are given by Ii3. 28\) . and the optimal stopping boundaries p^ 
and g*, such that 0<p* <p< l/2<g<g^, < 1, are uniquely determined by the coupled 
system of the equations in ^3.21^ with Daij)) and Dai^q) , i = 0, 1, given by \3. 25\) - l[3.26\) . where 
the functions Guiji) and Sii^n) , i = 0,1, are defined in l{3.15\) - l[3.18\) and ^3.2^^ , respectively. 

Proof. Since the verification of the assertions stated above can be done using similar ways 
of arguments, we present the proof of the second one only. Namely, we show that the functions 
in f l4.3p and f l4.4p coincide with the value functions in fl2.18p and fl2.19p . respectively, and the 
stopping times Q , i = 0, 1, from fl2.20p and thus the switching times (C*„)ngN from f l2.2ip are 
optimal with the boundaries p^. and g* specified above. For this, let us denote by Uq{tx) and 
Ui{'k) the right-hand sides of the expressions in fl4.3p and (14. 4p . respectively. Hence, applying 
Ito's formula to e^'''*f/j(nf ■*) , i = 0,1, and taking into account the smooth-fit conditions in 
flQij) . we obtain: 

e-^'-* U,{UP) = U,in) + r {UU, - r[/,)(n«) /(n« ^ p„ n« ^ g,) ds + M« (4.5) 

^0 

where the processes M(*) = (Mf ^)i>o defined by: 

Mf^ = / e-'' U[{Iif) n«(l - n«) dBs (4.6) 

Jo 

are continuous square integrable martingales under the probability measure with respect 
to the filtration (J-t)t>o, for every i = 0, 1. The latter fact can easily be observed, since the 
derivatives U-^n), i = 0,1, are bounded functions. 
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Taking into account the assumptions in fl3.32p . it is shown by means straightforward com- 
putations and using the properties of the functions t/i(vr), i = 0,1, that the conditions of 
fl2.36p and (12.371) hold with < < p and g < g* < 1. These facts together with the 
conditions in fl02D - flZ^ and fl^:^ yield that the inequalities (LqUo - rUo)in) > -(1 - vr) 
and (Lif/i — rf/i)(7r) > — tt hold for all vr G [0, 1] such that tt 7^ and tt ^ q^,, as well as 
^o(^) ^ an -\- 1/1(71) and f/i(7r) < 6(1 — vr) + f/o(7r) are satisfied for all tt G [0, 1]. It also follows 
from the regularity of the diffusion processes 11'^*^ i = 0,1, solving the stochastic differential 
equations in (12 .Mp and (12.150 . that the indicator which appears in the formula (14.50 can be 
ignored. We therefore obtain from the expression in (14. 5 p that the inequalities: 

-•Co 



a e 



-rco n(o) + f ° e-"^ (1 - ) ds + e-'"^" UiiufJ) (4.7) 
"Jo ° 



> e-'-^" UoiufJ) + e-'-^ (1 - nf ) ds > Uoin) + 
Jo 



(0) 
Co 



be-"-^' (1 - n[f ) + [' e-'' ds + e-'^' Uo{Iif) 
^ Jo 

>e-"^^f/i(n[f)+ / e-''^n«cis > f/i(7r) + M^'^ 
^ Jo ^ 



(4.8) 



hold for any stopping times Q of the processes 11*^*^ , i = 0,1, respectively. 

For every i = 0, 1, let (xj „)„gN be an arbitrary localizing sequence of stopping times for 
the processes M^'^^ . Then, taking the expectations with respect to the probability measure Pn- 
in (I4.7p -( l4l8l) . by means of the optional sampling theorem (see, e.g. [211 Theorem 3.6] or [20t 
Chapter I, Theorem 3.22]), we get: 



En 

>En 



CoA>co,7i 







(4.9) 



CoA^co.i 



e-'^"(l-nl°hc/s 



l(i) 

ClA><i,. 



ClAxi,, 



(4.10) 



ClA>fl,n 







for all TT G [0,1]. Thus, letting n go to infinity and using Fatou's lemma, we obtain that the 
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inequalities: 



En 



a e 



Co 



be-"^^' (1 - n[f ) + / e-'' ds + e~'<' Uoiuf^) 
Jo ^ 



> f/o(vr) 
>U,{tt) 



(4.11) 
(4.12) 



are satisfied for any stopping times Q, i = 0,1, and all vr G [0, 1]. By virtue of tlie structure 
of tlie stopping times in fl2.20p . it is readily seen that the equalities in (14.111) and (I4.12p hold 
with Q instead of Q, i = 0,1, when either vr < p^. or vr > , respectively. 

It remains to show that the equalities are attained in (14. lip and (14.12^ when (* replaces 
(i, i = 0, 1, for < vr < g*. By virtue of the fact that the functions Ui{n) , i = 0,1, with the 
boundaries and g* satisfy the conditions in (I2.32p and (I2.33p . it follows from the expression 
in (14.51) and the structure of the stopping times in (I2.20p that the equalities: 



:i-nl°))rfs = f/o(7r) + M;° 







(0) 

Co A>fo,„ 



Ci Aj<i,„ 



e-^^n«c?. = ^i(7r) + M«^,„ 



(4.13) 
(4.14) 



are satisfied for all tt G [0,1]. Observe that the integrals here are finite (Pj^-a.s.) as well as 
the processes (M^!^^J(>o, i = 0,1, are uniformly integrable martingales. Therefore, taking 
the expectations in (I4.13P and (14.140 and letting n go to infinity, we can apply the Lebesgue 
dominated convergence theorem to obtain the equalities: 



En 

En 



Co* 



^.rs _ n(o)) ^5 + e-r-cs UiiufJ] 



'J 



f/o(vr) 
f/i(7r) 



(4.15) 
(4.16) 



for all 7r G [0, 1]. The latter, together with the inequalities in (14. lip and (14.121) . directly imply 
the desired assertion. □ 

Remark 4.3. The results formulated above show that the following sequential procedure is 
optimal. Being based on the observations of X, we construct the posterior probability process 
n and stop the latter for the first time as soon as it exits either the region {^g^, h^) or (p^,, g^,) , 
appropriately, and then conclude that the process has switched either from to 1 or from 
1 to 0, respectively. Then, we continue to observe the process 11 which is currently located 
either in the regions [0, and [/;,*, 1] or in [0,p*] and [g*, 1] , and stop the observations as soon 
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as it comes to the opposite region. We may thus conclude that 6 should have switched either 
from to 1 or from 1 to 0, respectively, and continue the procedure from the beginning. 

Remark 4.4. Taking into account the results obtained above, we may also conclude that 
the appropriate minimal Bayesian risk functions take the form: 

V*{7r) = mm{Vo{7T), V^*{7v)} and U*{n) = min{f/*(7r), U^in)} (4.17) 

for TT G [0,1], where the functions V*{n) and U*{n), i = 0,1, are defined in fl2.8p -( l2l9|) and 
fl^:Tg]) -f iro]) . respectively. It is also seen that if either V*{7r) = V*{tt) or U*{tt) = U*{tt) holds 
for any tt G [0, 1] fixed, then the sequences (T*„)neN or given by f l2.13p and fl2.2ip are 

optimal in (I4.17p . appropriately. 
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